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0> ; ABSTRACT 

^5 , We present a measure of the inclination of the velocity ellipsoid at 1 kpc below the 

Galactic plane using a sample of red clump giants from the RAVE DR2 release. We 
find that the velocity ellipsoid is tilted towards the Galactic plane with an inclination 
of 7.3 ±1.8°. We compare this value to computed inclinations for two mass models of 
■ the Milky Way. We find that our measurement is consistent with a short scale length of 

r> \ the stellar disc [Rd — 2 kpc) if the dark halo is oblate or with a long scale length (Rd — 

^ . 3 kpc) if the dark halo is prolate. Once combined with independent constraints on the 

flattening of the halo, our measurement suggests that the scale length is approximately 
halfway between these two extreme values , with a preferred range [2.5-2.7] kpc for 
a nearly spherical halo. Nevertheless, no model can be clearly ruled out. With the 
continuation of the RAVE survey, it will be possible to provide a strong constraint on 
the mass distribution of the Milky Way using refined measurements of the orientation 
of the velocity ellipsoid. 

Key words: Stars: kinematics ~ Galaxy: fundamental parameters - Galaxy: kine- 
matics and dynamics. 



1 INTRODUCTION 

Our understanding of Galactic stellar populations and kine- 
matics makes regular progress with the advent of new large 
Galactic stellar surveys providing distances, photometry, 
radial velocities or proper motions. Our Galaxy is at the 
present the only place where we can probe the 6D phase 



space of stellar positions and velocities. For instance, the 
Galactic 3D potential c an be probed through the orbits of 



the Sagittarius stream (llbata et al. 2001 : Read & Moore 


2005 


; iNewberE et al. 1 12OO2I: Fellhauer et al. 1 20061: 


Helmi 


2004 


) or Palomar 5 tidal tails ([Odcnkirchen et al. 


' '2003|; 


Grillmair & Dionatos 1 l2006l: iGrillmair & Johnson 


200i), 



or th rough the kinematics of halo stars l|Battaglia et al. 1 
I2OO5I ). At smaller scales, the potential can also be anal- 
ysed through the force perpendicular to the galactic 
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plane (Oort 1960; Creze et al. 1998; Kuijken & Gilmore 
1989a,b,c,1991; Siebert et al. 2003; Holmberg & Flynn 
2004) or through the coupling between the 3 component s 
of the velocity in the solar neighbourhood (|Bienavme II 19991 ). 

Here, we concentrate on the question of the orienta- 
tion of the velocity ellipsoid that is known to be tightly 
relat ed to the shape and symmetry of the galactic poten- 
tial (lOllongren Ill962l:[ Hori &: Lui 1 1 19631 : iLvnden-Bellll 19621 : 
lAmendt fc Cuddeford Ill99ll ). 

In spite of the long interest in this problem, measuring 
observationally the orientation of the velocity ellipsoid out- 
side of the galactic plane has proven to be very difficult. This 
is due mainly to the absence of reliable distances away from 
the Solar neighbourhood. Despite this limitation, the first 
stellar stream detected within the Milky Way halo towards 
the n orth Galactic pole by iMaiewski. Munn fc Hawlev I 
ljl996h shows a velocity tilt, the ellipsoid being inclined to- 
wards the Galactic plane. This tilt could result from the 
expected velocity correlation induced by a spher oidal poten- 
tial if these stars had similar integrals of motion (jBienavme I 
1 19981 ). However we note t hat this stream is no t detected 
locally in the RAVE data (|Seabroke et al. 1120081 ). 

Building realistic Galactic potentials shows that the 
main axis of the velocity ellipsoid, at 1 kpc above the 
Galactic plane, points in the direction of the z-axis of 
symmetry of the Galaxy towards a point located at 5 to 8 
kpc behind the Galact ic centre : for in st ance from numerical 
orbit computations l|BinnevJ 19831: Kuiiken fc Gilmore I 
Il989al ) or applying to the ICarlberg fc Innanen (|l987l ) 
Galactic potential the lAmendt fc Cuddeford (1991) formu- 



lae. Such estimates of the velocity ellipsoid tilt are necessary 
for an accurate determin ation of the asymmetric drifu 
l|Binnev fc Tremaine |[l987h . and for a correct measurement 
of th e force perpendicular to the Galactic plane (|Statler I 
1 19891 ). 

In this paper, we study the 2D velocity distribution 
perpendicular to the Galactic plan e for a sample of red 
clump stars from th e RAVE survey jSteinmetz et al."ll2006l : 
IZwitter et al. |[2008l ). These stars are selected between 500 pc 
and 1500 pc below the Galactic plane and provide a mea- 
surement of the tilt of the velocity ellipsoid at ~1 kpc. In 
Section [2l we present the selection of the sample while Sec- 
tion |3] focuses on the measurement of the inclination and 
possible biases. Finally, in Section |4] we compare our mea- 
surement to computed inclinations for two extreme classes 
of mass models and we discuss possible outcomes of this 
measurement. 



2 SELECTION OF THE SAMPLE 

Our sample is drawn from the second data release of the 
RAVE survey (|Zwitter et al . 2008) containing about 50 000 
stellar radial velocities and 20 000 measurements of stellar 
parameters. We focus on red clump giants towards the South 
Galactic pole to maximize the distance from the plane and 



^ The asymmetric drift is the tendency of a population of stars 
to lag behind the local standard of rest for its rotational velocity, 
the lag increasing as a function of age. 



minimize the interstellar extinction. Hence, we select our 
targets in a cone with b < —60°, and we use a colour- 
magnitude criterion following IVeltz et al. I l|2008l ) to select 
our candidate red clump stars: 2MASS J — K colour within 
0.5-0.7 and K < 9.3. 

This colour-magnitude cut selects mainly red clump 
stars whose luminosity function (LF) is well defined and 
approximately Gaussian: Mk = —1.6 ± 0.03. Also, the red 
clump LF is narrow, the dispersion of the Gaussian LF 
being 0.22 mag in t he if-band, a nd is nearly independent 
of the metallicity lAlves I I2OOOI ). It makes this popula- 
tion particularly suited to study the kinematics of stars 
away from the solar neighbourhood as reliable distance 
estimates can be obtained. Also the extinction in the 
if-band remains low, <Ak>= 0.007 mag with a maximum 
extinction of Ak = 0.05 mag for t his region of the sky 



(Schlcgc l. Finkbeiner fc Davis I Il998l ). Hence, extinction 
does not contribute significantly to our error budget: the 
average error on the distance is less than 1% with a maxi- 
mum value of ~2% for the limiting magnitude of our sample. 

The selection criterion, retaining only the objects with 
a proper motion value in the RAVE catalogue, restricts the 
sample to 763 red clump candidates spanning a distance 
interval from the Sun of 500 to 1500 pc. A small fraction 
of these selected stars are dwarfs or subgiants. According 
to the photo metric and kinema tic modelling of the SGP 
and NGP bv lVeltz et al. ] l|2008l ). we can estimate that, at 
the limiting magnitude of our sample mK=9.3, 75% of the 
sample are red clump stars, 10% dwarfs and 15% subgiants. 
Brighter than this limit, the fraction of red clump stars 
is larger and the quoted fractions are upper limits to our 
contamination fraction. 

We clean our sample further using a kinematic se- 
lection. We select stars based on their velocities with the 
following criteria VU^ -I- -|- < 200 km s"^ and 
V < 100 km s~^. This selection enables us to remove 
the nearby dwarfs whose distance is overestimated by a 
factor 14 due to their fainter absolute magnitude, hence 
an overestimation of their velocities. The resulting sample 
contains 580 red clump candidates in the direction of the 
South Galactic pole whose distribution in velocity space 
is presented in Fig. [1] In this figure, the contours depict 
the distribution of the original sample smoothed by the 
individual errors while the dots show the location in velocity 
space of the remaining 580 stars after the velocity selection. 

We test our selection criteria using the second year 
observation from the RAVE survey (as a reminder, RAVE 
DR2 contains the first year of observation -i.e. DRl- and 
the second year of observation). For these objects, RAVE 
provides measurements of the stellar parameters including 
an estimate of the gravity. The sample selected from second 
year data contains 294 stars with logg measurements, with 
231 stars matching the velocity criteria. The histograms of 
log g for each subsample are presented in Fig. [2] where the 
black histogram presents the distribution of log g for the 294 
second year stars and the dashed histogram the subsample 
matching our velocity criteria. The red clump giants span 
a large range in gravity depending on their metallicity: 
logg = 2.08 for the metal-poor, low mass end and reaches 
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Figure 1. Selection of the red clump sample in U , V, W velocity space. The contours show the distribution of the 763 red clump 
candidates belonging to the original sample, smoothed by the individual errors, while the dots represent the location in the velocity space 
of the 580 stars in the final sample. The contours encompass 90, 70 , 50 and 30% of the total sample. 
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Figure 2. Top panel: histogram logg for the 294 stars with logg 
measurements in RAVE DR2. Black line full subsample, dashed 
line subsample matching the velocity criteria. The red clump gi- 
ants cover the region in log g = 2 — 3 depending on the mctallicity 
or mass. A conservative estimate of RAVE standard error on logg 
is 0.5 dex. Bottom panel: fraction of stars rejected by the velocity 
criterion as a function of logg. 



up to lo gg = 3 for the high mass, metal-rich red clump 
objects (|Zhao. Qiu fc Zhang I |2000| ). This figure clearly 
indicates that our velocity criteria is efficient for rejecting 
dwarf stars (with high logg) but also removes a small 
fraction of stars with lower logg, primarily subgiants and 
giants on the ascending branch and also a few red clump 
stars. Nevertheless, these objects have large velocities and 
fall in the tails of the velocity distribution. Therefore, they 
affect only marginally the measurement of the inclination, 
as our measurement is driven by the larger number of stars 
in the bulk of the velocity distribution. 

It is worth noting that due to the uncertainties in RAVE 
logg measurements whic h are 0.5 dex for a typical RAVE 
star (jZwitter et al. |[2008l ), it is not possible to obtain a firm 
estimate of the contamination in our sample, nor to use the 



RAVE log g estimates to refine our sample. Also, log g mea- 
surements are only available for less than half of our sample 
as stellar parameters can not be estimated from the spectra 
collected during the first year of operation of RAVE. Nev- 
ertheless, considering a 0.5 dex error on logg, we estimate 
that the contamination using the velocity criteria is reduced 
to ~10% which is to be compared to more than 20% without 
the velocity criteria. We will detail the effect of this contam- 
ination on our measurement in the next section. 



3 MEASURING THE TILT 

The tilt angle S of the 2D velocity distribution is given by 
the relation: 



tan 25 — 



(1) 



where amy, Ou and (jy^ are the velocity distribution 
moments. In the local velocity coordinates, the velocity in 
the radial direction is given by the U component of the ve- 
locity vector (positive towards the Galactic centre) and the 
vertical velocity by the W component of the vector positive 
towards the North Galactic pole, while the V component 
is positive towards the Galactic rotation (not used in Eq.[T|. 

The computation of the inclination is straightforward 
for a sample with small and homogeneous errors. Neverthe- 
less, to lower the effect of foreground dwarfs and giants and 
the contamination due to high-velocity stars, we make use 
of a velocity cut-off to select our sample. In this case, as our 
errors in the U and V velocity directions are large, a direct 
measurement of the tilt angle may be subject to bias and our 
selection criteria must be studied as our error budget may 
not be dominated by the size of the sample (see Section rS.lf) . 

Also, the local velocity ellipsoid is not a smooth dis- 
tribution and clumps are present on both s mall and large 
scale s in the velocity space (see for e xamp le 
19981: IChereul. Creze fc Bienavmil Il998l : 



Dehnen II2OO0I. 



Famaev et al. I 



2005h . These substructures prevent determining the age- 
velocity dispersion relation in the U and V directions 
ISeabroke fc Gilmore I |2007| ) and may also influence the 
measured tilt angle. We will discuss the effect of such sub- 
structures on our measurement in Section [3.21 
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Figure 3. Distribution of the velocity errors in our sample; full 
line U component, dashed line V component, dotted line W com- 
ponent. The difference between this three distributions arises due 
to the relative contribution of proper motion and distance errors 
to the radial velocity errors for the velocities along the three car- 
dinal directions. 

Finally, if our selection criterion is efficient at rejecting 
the foreground stars, only the tails of the velocity distribu- 
tion are affected by the velocity cut-off. As their space veloc- 
ities are overestimated, such foreground objects will impact 
on the measurement of the inclination. We will discuss this 
particular point in Section [3.31 



3.1 The effect of errors and velocity cut— off 

Our velocity errors in the cardinal directions are not ho- 
mogeneous because the U and V components of the veloc- 
ity vector are dominated by the proper motion contribu- 
tion while the W component is primarily measured from 
the RAVE radial velocity. Fig. |3] presents the distribution 
of errors for our sample in the U , V and W components 
as full, dashed and dotted lines. It is clear that the mode 
of the velocity error distributions for the U and W com- 
ponents, the ones we are primarly interested in, differ by a 
factor 4: ~ 5 km s~^ for the W component while for the U 
component the distribution peaks at ~ 20 km s~^. 

This large difference results in an anisotropic smooth- 
ing of the observed velocity ellipsoid which, combined with 
our velocity criterion, biases the measurement of the in- 
cMnation towards a lower value. This bias is due to the 
structure of Eq. [T] where the error anisotropy results in 
an extra component E on measured velocity dispersions. 
If we consider only the extra term on the U component we 
have au (measured)^ — au{true)^ + and the cross-term 
auw (measured) = puwo'w \/o-u{true)'^ + E'^, puw being 
the correlation coefficient. With this notation, it is clear that 
the contribution of the additional error term is larger for the 
denominator than it is for the numerator, hence producing 
an underestimate of the true inclination. For comparison, a 
linear fit would not be biased due to the asymmetry of the 
errors but unfortunately it is more sensitive to outliers. To 



overcome this problem, we use a Monte Carlo sampling of 
the velocity error distributions. We add a random velocity 
term to the V and W components, degrading the accuracy 
of the two velocity components, so that the resulting error 
distributions match the U velocity error distribution. For 
the U velocity, it is randomly drawn from its original error 
distribution. This procedure enables us to obtain isotropic 
error distributions for all three components, degrading the 
two best distributions to the level of the least accurate dis- 
tribution. The inclination is then computed using Eq.[T]after 
applying the velocity criterion. 

We tested this procedure on a simple velocity ellipsoid 
model using the RAVE error laws and standard velocity dis- 
persions for the Galactic old disc population, leaving aside 
the thick disc: au ~ 31 and aw ~ 17 km s^^ . The size of the 
sample was set to 1000 data points and we varied the incli- 
nation of the ellipsoid from 1 to 20° . The results of this test 
are shown in Fig. |3] where the direct measurement is pre- 
sented as a dotted line and the Monte Carlo determinations 
by the open circles with error bars for one random realisation 
of a velocity ellipsoid. The one to one relation between the 
original and recovered angles is sketched by the dashed line. 
Below 2-4° for the inclination, depending on the realisation 
of the ellipsoid, both methods predict the same inclination 
but above this threshold, the Monte Carlo sampling recov- 
ers the proper value of the angle. On the other hand, the 
direct measurement, applying Eq.[T] always underestimates 
the true angle with a bias rising with the tilt value. This 
test clearly indicates that the Monte Carlo sampling of the 
errors is best suited to measure the tilt of the velocity ellip- 
soid, while direct measurements using Eq. [l]are subject to 
strong bias in the case of heterogeneous error laws. We note 
also that the value of the bias depends strongly on the ran- 
dom sampling of the ellipsoid, with a bias varying between 
2 and 4° at 7° . This spread in the bias value becomes larger 
as the tilt value increases and indicates that even with a 
proper model to estimate the bias, it can hardly be used to 
correct the direct measurement. 

The procedure is then applied on the RAVE sample and 
the result is presented in Fig. [5] which shows the distribution 
of inclinations in degrees obtained by sampling the error 
distribution 25 000 times. The mean inclination measured 
is 7.3° with a standard deviation of 1.8°. If the anisotropy 
of the error distributions were not taken into account, the 
measured inclination would have been 6.3° or 1.0° too low. 
The 2D representation of the velocity ellipsoid inclination is 
shown in Figure |S] The colour-coding follows the density of 
stars per bin in the region of the {U,W) space, where the 
2D distribution of the {U,W) velocities has been convolved 
by the individual errors. The measured inclination and l-cr 
errors are presented as white lines (full line for the mean 
value and dotted lines for the errors). 

3.2 Effect of substructures 

Substructures such as the Hyades, Pleiades or the Hercules 
groups are well known features of the local velocity ellipsoid, 
and are easily seen in the velocity sp ace obtained from the 
Hipparcos mission (s ee for e xample Famaev et al. 20051: 
Dehnen fc Binnev I 1 19981 : IChereul. Creze fc Bienavme I 
19981 ). These structures may have a wide range of origins 



such as the disruption of clusters, resonances associated with 
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Figure 4. Test of the Monte Carlo method to measure the tilt 
angle of the velocity ellipsoid. The original versus recovered tilt 
angle are presented for one realisation of a velocity ellipsoid hav- 
ing (7(7 = 31 and criy = 17 km and sampled using 1000 data 
points. The RAVE error laws for the velocities are used. The 
dashed line shows the one to one relation while the dotted line 
is a direct measurement using Eq. [T] The results from the Monte 
Carlo sampling of the error laws are depicted by the open circles 
and the error bars are the standard deviation of 5000 resampling 
for each value of the tilt angle. 
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Figure 5. Distribution of the measured inclination of the velocity 
ellipsoid per 0.05° bin. This distribution is obtained using a Monte 
Carlo sampling of the error distribution. The mean inclination is 
found to be 7.3° with a standard deviation of 1.8°. The grey line 
is a Gaussian function with identical parameters. 



the bar or spiral arms fsee for examDle iFamaev et al. 1 20051: 


Famaev. Siebert & Jorisscn 1 l2008l: Minchcv & ( 


3uillcn ' 


2008; Dehncn 


2000: Dc Simone. Wu & Trcmainc 


2004). 



Nevertheless, the average velocity error in our sample does 
not allow us to distinguish these substructures. 

To test the influence of these velocity substructures 



The tilt of the velocity ellipsoid from RAVE 5 

200 



100 




-100 



-200^^^^^^^^^^^^^^^^^^^^™ 

-200 -100 100 200 

U (km/s) 

Figure 6. Velocity distribution in the {U,W) plane from our 
sample after sampling the error distribution. The measured incli- 
nation and l-CT range are presented by the full and dashed white 
line. The colour— coding follows the density per bin. 



Table 1. Fra c tion and number of stars in structures in the 
I Famaev et al. I l|2005l ) sample as a function of height above the 
Galactic plane. 



z Ntot Fraction in Num. in Num. in Num. in 

(pc) structures Hya/Plei Sirius Hercules 



0-100 1361 0.40 127 63 125 

100-200 1337 0.33 111 49 131 

200-300 844 0.28 46 27 94 

300-400 422 0.22 1 20 37 

400-500 177 0.19 7 8 



on the tilt de te rmina tion, we use the local sample from 
iFamaev et al~l l|2005h . This sample provides not only 
accurate velocity vectors for about 6500 stars in the solar 
neighbourhood, it also provides an estimate of the relation 
of a star to the identified velocity substructures. This allows 
us to separate the background ellipsoid from the known 
overdensities. 

W e first estimate t he fra ction of stars in substructures 
in the IFamaev et al. I |200^ sample as a function of z, 
the height above the Galactic plane. The number of stars 
in structures is larger closer to the plane: 36% of the 
stars are in structures in the 0-200 pc interval while 25% 
are found in structures between 200 and 500 pc. This 
drop in number of objects in structures is sharp as seen 
from Table [T] the fraction in objects in structures being 
lowered by over a factor 2 between and 500 pc. As 
our sample covers a distance below the plane from 500 
to 1500 pc, we extrapolate this behaviour at higher z to 
estimate the contamination arising from substructures 
in our sample. Using a conservative extrapolation, we es- 
timate the contamination in our sample to be lower than 7%. 

In a second step, we test the influence of the substruc- 
tures on the measur ed inclination. We can not use the 
IFamaev et al. I (|2005l ) sample directly, as the presence of 
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Figure 7. Influence of st ars in velocity groups on the measured 
inclination. A sample of iFamaev et al. I ||2005| ) stars belonging 
to groups is randomly added to the RAVE sample, varying the 
contamination from 1 to 20%. The tilt is measured following the 
same procedure as for the pure RAVE sample. The thick line 
represents the average deviation in degrees {Stat = ^measured ~ 
Strue) found for 25 000 repeats per contamination fraction. The 
dashed line is the standard deviation of the repeats. 



pertubations in the plane makes any attempt to disentangle 
the effect of groups from the effect of the pertubations on 
the inclination hazardous. Therefore, we proceed as follow. 
We add an additional population, drawn from a subset of 
the Famaey et al. sample, to our RAVE sample. This subset 
is randomly selected from the set of stars belonging to 
groups in the distance interval 300 to 500 pc. We add it to 
the RAVE sample varying its fraction relative to the RAVE 
sample from 1 to 20%. This procedure enables us to mimic 
as closely as possible the velocity distribution of the groups, 
which is not homogeneous and strongly varies as a function 
of distance to the plane. This operation is repeated 25 000 
times for each fraction of the contamination to ensure a 
proper coverage of the possible cases. 

The results are shown in Fig. [7] where the mean de- 
viation {Smeasured — Strue) lu dcgrces as a fuuctiou of the 
fraction of stars in groups is drawn as a thick line. The 
standard deviation of the repeats is shown as dashed lines. 
We note that the average tilt for the Famaey group sample 
used here is 0.07 ± 2.45°, very different from the value of 
the inclination in our sample. From this figure, we see that 
the influence of velocity structures on the measured tilt is 
low, around —0.03° for a contamination of 7% with a stan- 
dard deviation below 0.1°. 7% being an upper limit for the 
contamination in our sample, we do not expect groups to 
affect our measurement of the tilt. Indeed, the mean devia- 
tion combined with the standard deviation measured from 
this experiment contributes to not more than 0.1°, less than 
6% of our estimated errors. 



3.3 Effect of foreground stars 

If, as we saw above, the velocity structures do not influence 
largely the orientation of the velocity ellipsoid, the fore- 



ground stars (dwarfs, subgiants and giants on the ascending 
branch) can be more problematic. The fact that their 
velocities are overestimated in the U direction, because of 
the overstimate of their distances, will add a component 
with low inclination to the observed ellipsoid. The bias due 
to these objects is an underestimate of the tilt at a given 
distance. 

The contamination by foreground objects is about 10% 
in our sample (see Section and a factor 14 overestimation 
of the U velocity for the dwarfs will render their velocity 
ellipsoid almost uniform in the velocity interval we consider. 
The velocity dispersion of this foreground population will 
be large, over 400 km.s~^ instead of ~31 km.s"^ for au, 
due to the distance overestimate. For the other sources 
of contamination (subgiants and giants on the ascending 
branch), the distance is overestimated by a factor 2 or less, 
and their impact on the velocity ellipsoid is lower. 

To obtain an upper limit of the effect of the foreground 
population we rely on a resampling technique, replacing 
10% of the sample by a random realization of a thin 
disc population with no inclination. The overestimate of 
the distance is then translated into an overestimate of 
the velocities, and the final inclination of the velocity 
ellipsoid is measured applying the same procedure as 
above. The difference between the distribution with and 
without resampling provides an upper limit on the effect of 
foreground objects on our measurement. We note here that, 
if adding a population with no tilt is in principle similar 
to the experiment done in Section [3.21 here we incorporate 
the distance overestimation. Furthermore, the added test 
stars are not restricted to the region in velocity space of 
the groups. 

Figure [8] presents the results of this resampling. The 
black histogram shows the distribution of the 25 000 mea- 
surements while the grey Gaussian curve is the Gaussian 
representation of the distribution obtained in Fig. [5] The 
presence of a population with no inclination does produce 
an observable bias in this experiment: we observe a shift 
in the maximum of the distribution. Nevertheless, this bias 
is small, the measured offset is 0.15°, much lower than 
the standard deviation of the distribution while it is a 
worst-case scenario. Indeed, here we did consider only the 
contamination by foreground dwarf stars while our real 
contamination is a mixture of dwarfs and subgiants. In 
the latter case, the overestimate of the distance is much 
lower, as these stars have a mean distance from the plane 
that is larger. Hence, their impact on the velocity ellipsoid 
orientation will be lower than for dwarfs. We note also 
that the presence of a foreground population renders the 
distribution non-Gaussian, adding a tail to the low tilt 
angle part of the distribution which is not observed in Fig. [5] 

We can conclude that our estimate of the inclination of 
the velocity ellipsoid is robust and that the presence of a 
population of foreground objects does not introduce a sig- 
nificant bias in our measurement. Indeed, combining both 
the effect of foreground stars and of the possible velocity el- 
lipsoid substructures in a worst-case scenario, the resulting 
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Figure 8. Results from the resampling study. 10% of the red 
clump sample has been replaced by a population of foreground 
objects with no tilt and velocities in the U direction overesti- 
mated by a factor 14. The black histogram is the distribution 
of measured inclination for this new sample following the same 
procedure as for Fig. [5] while the grey Gaussian curve is the dis- 
tribution of inclination without the resampling from Fig. |5] 



bias amounts to ~10% of our errors. At this level, the biases 
do not affect our conchisions. 



4 RELATION TO THE MASS DISTRIBUTION 
IN THE GALAXY 

The tilt of the velocity ellipsoid is intimately linked to the 
mass distribution in the Milky Way and more specifically 
- if we trust our knowledge of the structure of the Galac- 
tic disc - to th e fjattening of the h alo. W e start from the 
mass model of iDehnen fc Binnev I l|l998lfl and its revised 
parameters provided by Binnev fc Tremaine I (|2008t ) in their 
table 2.3 (hereafter BT08). This revision proposes two mod- 
els, referred to as model I and II, which match both local 
an d non-local data. These m odels are modified versions of 
the lDehnen fc Binnevi (119981 ) models 1 and 4. 

As noted in these references, a crucial parameter is 
the scale length of the disc whose value lies in the range 
2-3 kpc. The two models are set on the upper and lower 
bound for this parameter. Model I presents a mass model 
with a short scale length [Rd = 2 kpc) that induces a 
strong contribution from the disc for the potential at the 
solar radius up to 11 kpc. On the other hand, model II 
has a larger scale length [Rd = 3.2 kpc) and therefore, 
the halo contribution to the rotation curve dominates 
at the Sun location and beyond. This is also seen from 
the global shape of the potential where for model II the 
isopotentials are more spherical than for model I (see 
figs. 2.19 and 2.21 of BT08). For a detailed description of 



^ The Galactic potentials are computed using the GalPot pro- 
gram written by W. Dehnen. This program is available within 
the NEMO package: [http: / /ca rma. astro. umd.edu/nemo/ 1 



Table 2. Modification to the mass models I and II of BT08 table 
2.3 used in Section |4] The halo density is given in M0pc-3. The 
modified models are built, modifying the halo parameters, to keep 
the rotation curve almost unchanged in the disc. The bulge and 
disc parameters arc fixed to the BT08 values. 



Phalo 



c/a Model I Model II 



0.6 


0.838 


0.327 


0.7 


0.765 


0.293 


0.8 


0.711 


0.266 


0.9 


0.670 


0.245 


1.0 


0.635 


0.229 


1.1 


0.608 


0.215 


1.2 


0.585 


0.204 


1.3 


0.566 


0.195 


1.4 


0.548 


0.186 


1.5 


0.534 


0.179 


1.6 


0.520 


0.172 


1.7 


0.510 


0.167 



these two models, the reader is referred to chapter 2 of BT08. 

We use both models to discuss below the implications 
of the tilt on the possible models for the mass distribution 
in the Milky Way, focusing on the flattening of the halo in 
the two extreme cases. We note here that the region above 
(below) the plane between 1 and 2 kpc is best suited to 
separate the two classes of models. Indeed, in this region, 
the variation of the angle between the Galactic plane 
and the normal to the isopotentials as a function of the 
minor-to-major axis ratio c/ap is maximum. Hence, we 
expect the difi^erence between the predicted tilt angle to 
be the largest in the same region. At larger distances from 
the plane, the potential becomes more spherical and the 
difference vanishes between the models in terms of variation 
of the potential and inclination of the velocity ellipsoid. 

To measure the tilt of the ellipsoid as a function of the 
halo flattening, we vary the density minor-to-major axis 
ratio c/op of the halo from 0.6 to 1.7 for each model, the 
halo density being described by the relation 

( m\-"h ( ^ m \ ('^h-0h) 

lo ( — I 

\ah/ 



p{R,z) 



Phalo 



1 + 



(2) 



where the flattening c/up enters the equation through the 
parameter m — R^ + /{c/ap)^ , R and z being the 
Galactocentric cylindrical coordinates, is a scale param- 
eter with p oc m~°'^ if m <C au and p oc m~^^ for large 
m. 

While changing c/op, we change the halo density in 
order to keep the rotation curve almost unchanged in the 
plane. The bulge and disc components are flxed to the values 
of table 2.3 of BT08, the solar Galactocentric radius= 8 kpc. 
Keeping the rotation curve and the disc/bulge parameters 
unchanged enables us to study the influence of the halo flat- 
tening on the shape of the velocity ellipsoid, as the contribu- 
tion to the radial force of each component remains largely 
the same for each case. The corresponding density of the 
halo for each mass model is reported in Table [S] 

The inclination of the velocity ellipsoid is computed for 
a given mass distribution using orbit integration. A sin- 
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gle orbit is integrated over 30 rotations using a 4th order 
Runge-Kutta algorithm. The initial conditions are drawn 
from a Shu distribution function matching the local data 
iBienavme |[l999 l . For each potential, the orbit library con- 
tains over 2 million orbits from which we randomly select 10 
points per orbit in the last 15 rotations. We further restict 
the orbit library to data points matching the interval in R 
and z of the RAVE sample, this reduces the size of our final 
libraries to 6.10'* to a few 10^ points per library. 

We measure the tilt using a Monte Carlo selection of the 
orbits, requiring that the distribution of the selected orbits 
matches the selection function of our RAVE sample in the 
{R,z) plane. Here R is the Galactocentric radius and z the 
distance above (below) the plane. This selection function is 
obtained by convolving the distribution of the RAVE sample 
in R and z by their errors for each star. This procedure 
ensures us that the spatial distribution in R and z of the 
RAVE sample is well matched by the orbit selection. We 
select 5 000 orbits using the spatial constraints, about 10 
times larger than the observed sample but about 10^ times 
less than orbit library size to minimize the probability of 
the same orbit to be selected twice, and the tilt is measured 
using the associated velocities and Eq. [1] The measurement 
is repeated 500 times to obtain the mean inclination and 
dispersion. This procedure is repeated for each orbit library. 
The convergence of this procedure is tested using 1 000, 
5 000 and 10 000 orbits points for the selection. The result 
shows a very good stability of the mean inclination, of the 
order of a few lO"'^ degrees, while the dispersion increases 
as the number of orbits becomes lower. For the model I 
with c/flp — 1.0, we obtain respectively an inclination of 
9.89 ± 1.09, 9.83 ± 0.44 and 9.85 ± 0.36 for 1 000, 5 000 
and 10 000 orbits, which indicates that the gain in precision 
above 5 000 orbits is limited as the computing time scales 
linearly with the number of orbits. 

The resulting measurements are presented in Fig. [5] left 
panel where the full horizontal line is our measurement from 
Section [3] and the horizontal dashed lines are the la limit. 
The two remaining curves correspond to the measurements 
obtained from our orbit analysis. The top curve is our 
prediction for the class of models I of BT08 and the bottom 
curve for the class of models II. For comparison, the 
right panel presents for the direct measurement, without 
correcting the tilt for the velocity error anisotropy. The 
circles and crosses correspond to Monte Carlo realisations 
of the RAVE sample using the orbit libraries where the 
RAVE velocity errors have been applied on the orbit library 
directly. Circles and crosses are respectively for the class of 
models I and II. The error bars are the standard deviation 
of 1 000 realisations. 

The classes of models show the same general gross 
properties. The predicted tilt rises as the flattening de- 
creases, reaching a maximum in the prolate halo region. 
This maximum is expected and varies depending on the 
details of each model. It is due to the fact that, when c/up 
becomes large, the potential becomes separable in cylin- 
drical coordinates. On the other hand, if c/ap approaches 
0, the problem reduces to the plane-parallel case and the 
expected inclination at 1 kpc is 5 ~ 3°. 

The two classes of models provide different estimates 



for the tilt of the velocity ellipsoid within the limits of our 
sample. The tilt variation versus c/ap is larger for prolate 
models than it is for oblate models. For a density flattening 
of 0.6, the difference is only 2°, while for slightly prolate 
models the difference reaches up to 5°. Also, in the oblate 
case, the expected inclination rises more quickly than for 
the prolate case. If one compares the predictions for the 
two classes of models to the measured inclination, a clear 
tendency is present: the more massive the halo is, the 
more prolate it must be to match the tilt of the velocity 
ellipsoid at 1 kpc. The two extreme cases "separate" around 
c/ap ~ 0.9: while for a massive disc c/ap ^ 0.9 is necessary 
to reproduce the tilt, c/ap ^ 0.9 is needed for a massive halo 
in the 1-a limit. More speciflcally, the measured orientation 
of the velocity ellipsoid is consistent with a short scale 
length of the disc if the halo is oblate, while in the other case 
- a scale length of the disc of the order of 3 kpc - the mea- 
sured value for the tilt implies that the halo must be prolate. 

Using a direct measurement, applying the velocity er- 
rors on Monte Carlo realisation of the RAVE sample has the 
benefit of reducing the errors by a factor \/2 (Fig. |9] right 
panel). Nevertheless, the bias increases with the inclination, 
see Fig. |4j which results in the difference between the two 
models being lower. This direct measurement indicates that 
low values of c/ap are marginaly inconsistent with the mea- 
sured tilt, with c/ap > 0.7 being prefered even so the same 
general conclusions hold for both modeling technique. We 
note however that applying the errors on the orbit library 
(direct method) and correcting the anisotropy of the veloc- 
ity errors (unbiased measurement) produces slightly differ- 
ent predictions for the tilt. This is partly due to the fact 
that the U and V velocities are computed from the knowl- 
edge of distances and proper motions. Hence, U and V errors 
increases with distance which is not taken into account in 
the Monte Carlo simulation for the direct measurement, the 
RAVE sample being too small to estimate properly the error 
laws as a function of distances. This effect is also reduced in 
the unbiased measurement but is still present since the cor- 
recting term is added to the true error, hence distant stars 
will still have on average large errors while nearer objects 
will have on average smaller errors. 

Further constraints on the minor-to-major axis ra- 
tio are also available from independent studies. For ex- 
the flattening of the dark halo has been esti- 
from the shape of the Sagittarius dwarf tidal 
A value c/ap > 0.7, with a pr eferred flat- 
c/ap 



ample, 
mated 
stream 
tening 
(2001, 
(.200a 



1, is obtained bv | Ibata et al. 



of 

and iMaiewski. Skrutskie. Weinberg fc Qstheimer 
using respectively carbon stars and M-giants 
from the 2MASS survey a l ong t he orbit of Sagittarius. 
Ijohnston. Law fc Maiewski I ()2005h gives even stronger con- 
straint 0.75 < c/ap < 1.1 at a 3-a level, with oblate haloes 
strongly fav oured i f prece ssion of Srg's orbit is considered. 
In contrast, iHelmi I ()2004l ) and lLaw. Johnston fc Maiewski! 
(|2005h demonstrated that only Galactic potentials with pro- 
late halos could reproduce the velocity trends in the lead- 
ing debris with a preferred axis rat io c/ap = 5/3. However, 
iLaw. Johnston fc Maiewski I (|2005l ) explored a wide variety 
of Galactic potentials but failed to find a single orbit that 
can fit both the velocity trends and the sense of precession. 
Looking back at the solar neighbourhood, if the Sgr 
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c/a c/a 

Figure 9. Left panel: inclination of the velocity ellipsoid as a function of the halo flattening c/op in the RAVE selection function. The 
full and thin dashed-dotted horizontal lines correspond to our measurement and error bars using isotropic error laws. The two dashed 
curves correspond to the class of model I (top) and II (bottom) in BT08 for which we varied the halo flattening. Right panel: same as 
left panel but for a direct measurement. The horizontal line is the direct measurement of the tilt without correction for the velocity error 
anisotropy. The circles and crosses are for Monte Carlo realisation of the RAVE sample using the RAVE velocity error laws for the class 
of model I (circles) and II (crosses). The error bars are the standard deviation obtained from 1 000 realisations. 



stream is orbiting i n oblate and spherical po t ential s, Law 
et al. (2005) and iMartinez-Delgado et all l|2007l ) both 
predict that the Sun is currently bathing in a stream of 
debris from Sgr, passing both inside and outside the solar 
circle. Models orbiting in prolate potential s are on the other 
hand i nconsistent with this predi ction. iBelokurov et al. I 
l|2006t ). iNewberg et all (|2006l ) and ISeabroke et al I (|2008| ) 
all provide strong evidence for the absence of Sgr debris 
in the solar neighbourhood. iFellhauer et al. I ({2006 ) argue 
that the origin of the bifurcation in the Sgr stream is 
only possible if the halo is close to spherical, as the 
angular difference between the branches is a measure of 
the precession of the orbital plane. This suggests that 
the absence of the Sgr stream near the Sun is consistent 
with nearly spherical and prolate Galactic potentials and 

seemingly inconsisten t with oblate po tentials. However, 

recently iRuzic ka. Palo us fc Theis I (|2007l ) studied the Mag- 
ellanic System - Milky Way interaction using test particle 
simulations and compared them to HI observations. They 
concluded that c/up < 1 values (oblate halo) are prefered 
and allow a better match to HI observations. 



In Fig. |9l the preferred region by most studies, 
0.75 < c/op < 1, does not permit us to set strong con- 
straints either on the flattening nor on the mass of the 
disc. The measured value of the tilt falls between the two 
classes of models in the allowed region, and the error bars 
on the RAVE measurement do not permit us to tighten the 
paramete r space reliably . For strongly prolate halos, as sug- 
gested by iHelmil (|2004l ). the class of models II is preferred 
while short disc scale length are marginally rejected at the 
2- g level. If one adopts the axis ratio c/flp = 1 as prefer red 



by Maiewski. Skrutskie. Weinberg fc Qstheimer I (|2003l ) or 
llbata et al. 1 l|200ll ). the value of the tilt is better recovered 
with a model whose scale length of the disc lies in the range 



Rd = [2.5 — 2.7] kpc. Nevertheless, at the 1-cr level, large 
and short values for Rd are permitted with this analysis. 

Various studies in the literature have used star counts 
to constrain the scale length of the thin disc. For example 
recently Juric et al. (2008) measured the scale length of 
the stellar disc and found Rd = 2.6 kpc (±20%) using 
SDSS data. Similarly, usin g data from the Bolo gna open 
cluster survey (BOCCE), ICignoni et al~l ll2008i) f o und a 
scale length in the range 2.25-3 kpc and Oiha I (|200ll ) 
found 2.8 kpc using the 2MASS survey. These values are 
in good agreement with our finding but have similarly 
large error bars. Refining our measurement will provide an 
independent constraint on the scale length of the stellar disc. 



5 CONCLUSIONS 

We measured the tilt of the velocity ellipsoid at ~1 kpc be- 
low the Galactic plane using a sample of red clump giants 
from the RAVE DR2 catalogue. We find its inclination to be 
7.3±1.8°. Estimates of the effect of contamination by fore- 
ground stars and substructures have been shown to be small 
and their effect on our measured value can be neglected. 

We compared this value to predictions from two ex- 
treme cases of mass models for the Milky Way proposed by 
BT08. In the case of a massive disc with a small scale length 
{Rd = 2 kpc), the inclination is compatible with an oblate 
halo whose minor-to-major axis ratio c/ap is lower than 0.9 
at the 1-a level. On the other hand, in the case of a mas- 
sive halo with large disc scale length {Rd — 3 kpc), prolate 
haloes are prefered with c/op ^ 1. When a direct measure- 
ment is used, low values for c/up can be marginally rejected, 
indicating that c/up > 0.7. 

When further independent constraints from previous 
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studies are considered, we find that an intermediate value 
for the disc scale length Rd ~ [2.5 — 2.7] kpc is preferred 
for a nearly spherical halo, but no extreme model can be 
clearly ruled out, due to our large error bars. This range 
is in good agreement with other studies relying on star 
count analysis and deep photometric surveys. Nevertheless 
these results have large error bars of the same order as our 
measurement and cannot be used to further constrain the 
mass distribution. 

RAVE continues to acquire spectra and this work relies 
on the second data release of the survey. So far RAVE has 
collected more than 200 000 spectra, 4 times the size of 
the sample used here. With the current observing rate, we 
can expect to multiply by ten the size of our sample in the 
coming years which will allow us to significantly reduce our 
error bars. By the end of the survey, we will be able to 
provide a new mass model for the Milky Way galaxy with 
a constrained scale length of the disc and minor-to-major 
axis ratio of the dark halo. 
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